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We measure the sweep-to-sweep autocorrelations of blocked loops below and above the deconfinement transition 
for SU(3) on a 16 4 lattice using 20000-140000 Monte-Carlo updating sweeps. A divergence of the autocorrelation 
time toward the critical /3 is seen at high blocking levels. The peak is near f3 = 6.33 where we observe 440 ± 210 for 
the autocorrelation time of 1 x 1 Wilson loop on 2 4 blocked lattice. The mixing of 7 Brown- Woch overrelaxation 
steps followed by one pseudo-heat-bath step appears optimal to reduce the autocorrelation time below the critical 
/3. Above the critical (5, however, no clear difference between these two algorithms can be seen and the system 
decorrelates rather fast. 



1. Introduction 

In the Monte Carlo evaluation of physical ob- 
servables, we should gather independent samples 
to get statistically meaningful results. We gen- 
erate a sequence of configurations and then we 
choose samples separated by some sweeps to en- 
sure their independence. This separation should 
be larger than the Monte Carlo autocorrelation 
time t. Otherwise statistical errors should be 
corrected to take correlations into account. It is 
crucial in lattice QCD to estimate correctly these 
errors and the underlying autocorrelation time. 



In this paper we report autocorrelation measure- 
ments in pure SU(3) lattice gauge theory. [[[) 
The aims of our study are 

1. to extract autocorrelation time 

2. to see operator dependence of autocorrela- 
tion,!. e. blocked loops 

3. to see (3 dependence of autocorrelation 
which tells us some information about the 
critical slowing down toward the continuum 
limit. 
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Fig.l MC history of Wll at level=2-4, /3=6.38. 

2. Overrelaxation, Blocking scheme and 
Numerical simulation 

To update the gauge field here we use Cabbibo- 
Marinarig pseudo-heat-bath(PHB) with 3 SU(2) 
subgroups and to perform the overrelaxation(OR) 
we adopt Brown- WochQ microcanonical method 
with 2 randomly chosen SU(2) subgroups. 

We introduce a parameter K which denotes the 
mixing ratio between OR and PHB, namely we 
choose OR or PHB stochastically with the prob- 
abilities K/(K + 1) and 1/(K + 1). 

For blocking we use the SwendsenQ factor 2 
blocking scheme. For details and results of our 
MCRG study, see Ref.[| 

We work on a 16 4 lattice and 20000-140000 
sweeps are performed at all f3. All data are K = 9 
case except for (3—6.3 and 6.8. At /3=6.3 and 6.8, 
K=0, 3, 5, 7, 9, 12 were studied to check the 
efficiency of OR algorithm. 

After a few thousands of sweeps for thermal- 
ization, we measure operators: Wll, W12, W22, 
Wchair, Wsofa and Wtwist at every blocking level 
from 1 to 4 (i.e. from 16 4 to 2 4 ). Blocking is per- 
formed every 10 sweeps. All our simulations have 
been performed on 64-node AP1000 at Fujitsu. 

3. Results 

In Fig.l, we show typical Monte Carlo history 
of Wll (lxl wilson loop) at blocking levels 2-4 at 
0=6.38. 
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Fig. 2 p(t) vs. t, same data as in Fig.l. 

In this case no clear signal of sweep-to-sweep 
correlation at low blocking level can be seen 
whereas at high blocking level some structure 
emerges which we cannot see at lower levels. 

Then we define the autocorrelation function 
p(t) of observable O as 

(t)= ((o(i)-QA)m+t)-o B )} 

^(( (i)-o A )i)((o(i + t)-o B y) 

where ( ) denotes average over i and Oa = (0(i)) 
and O b = (0(i + t)). 

Fig. 2 shows the autocorrelation function p(t) 
for t=0-1000 sweeps for the same data as in Fig.l. 
In Fig. 2 we see long range autocorrelation, as 
stated above, at high blocking level. 

This tendency is very natural. Since Wll at 
high blocking level has effectively large area, then 
more Monte Carlo sweeps are needed to alter its 
value. In other words, Wll at high blocking level 
is not easily influenced by short range fluctuations 
whereas at lower level it is. 

Now we introduce integrated autocorrelation 
time Ti n t defined as: 

N N-t 
r int = p (0) + 2j2p(t)^ r (2) 
t=i * 

In this analysis we choose iV=1000, so that Ti nt 
can extract long range information. For the error 
estimation of correlated samples, Ti nt is the rele- 
vant quantity and enters in our final error analy- 
sis. 
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Fig.3(a) p{t) vs. t, for if =0,3,7,9 at /3=6.3. 

We have done this analysis not only for Wll 
but also for other measured loops such as W12, 
W22, Wsofa etc. There exist, of course, differ- 
ences among them about autocorrelation func- 
tion and time, however, we found that for a given 
blocking level, these differences are not signifi- 
cant. So from now on, we concentrate on Wll at 
level=4. 

As stated above, we have checked the efficiency 
of OR algorithm changing the mixing parameter 
K. 

In Fig.3(a), p(t) for t=0-600 at /3=6.3 for 
.K=0,3,7 and 9 are displayed. In this figure, K=7 
seems to work well to reduce the autocorrelation. 
To see this quantitatively, we show Tint defined by 
eq.(2) for several values of if at /3=6.3 in Fig. 3(b). 
We can see the efficiency of K=7 case also in this 
figure. 

From this analysis, we can extract the typical 
autocorrelation time of Wll at level=4 at f3—6.3 
such as: about 300 sweeps for PHB, about 100 
sweeps for OR at K=7. In ref. |J Decker and de 
Forcrand give similar numbers. These numbers 
are somewhat larger than those obtained by other 
authors. This is because we are using blocked 
loops as observables. To see the long range au- 
tocorrelation, extended operators look suitable 
because small operators are easily influenced by 
fluctuations. 

When we want to include error bars, we must 
be careful to control these long range autocorre- 



Fig.3(b) Ti n t for several K values at /3=6.3. 

lation effects. 

Let Spo(t) be the error on p(t) calculated as 
if samples were independent. Then if we con- 
sider the dependence between samples, the error 
is modified as 



(5 P {t)f = (5po(i)) a 



(3) 



This error on p(t) generates an error on the 
extracted autocorrelation times. We estimate it 
by 



_ Sp(T int ) 

"Tint — — 7 T T ii 

P{Tint) 



(4) 



To see the definition dependence of error anal- 
ysis, we also evaluate errors using jackknifc 
method. To perform jackknifc analysis for error 
estimates of autocorrelation function, we first di- 
vide the sample into bins. We get Sp(t) and then 
we can estimate the error of autocorrelation time 
following eq.(4). It is thought that if we dou- 
ble the size of a bin, error increases roughly by 
a factor \[2 when size of a bin is smaller than 
the autocorrelation time and then levels off when 
size of a bin exceeds the autocorrelation time. In 
our analysis there are many cases when the er- 
ror won't level off. In such cases the error on the 
error of autocorrelation function is very large. 

We show the results of this jackknife study also 
in Fig.3(b). The errorbar with full circle comes 
from using eq.(3) and that with open circle comes 
from jackknife analysis. In all cases of K, jack- 
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Fig. 4(a) p(t) vs. t, for several (3 at K=9. 

knife method estimates smaller errorbars than 
eq.(3). 

If we include the error and consider again the 
efficiency of OR algorithm, we can only say that 
K=7 looks optimal. To adopt OR algorithm, it 
is important to choose the optimal value of K 
to update the system efficiently. The optimum 
range of K is surprisingly narrow. 

We have performed the same analysis also at 
/?=6.8. At this (3, Tint is rather small for any 
K (typically less than 50 sweeps) and no clear 
advantage of OR can be seen. 

Finally we consider the (3 dependence of p(t) 
and T in t ■ 

We show the behaviors of p(t) for different val- 
ues of (3 in Fig. 4(a). We can see long range auto- 
correlations at ,3=6.33 and 6.43, whereas rather 
short range ones at /3=6.10 and 6.80. 

The autocorrelation time as a function of (3 
is displayed in Fig. 4(b). Tint seems to diverge 
around some (3, associated with the deconfine- 
mcnt transition. This indicates that around this 
(3 T we must be careful to evaluate the statistical 
average of physical observables as we stressed in 
the beginning of this paper. 

From these results, it is clear that if we want 
to stay in the confinement regime, we must stay 
below (3 T = 6.34 ± 0.03 on a 16 4 lattice. 

Combining with renormalization group analy- 
sis at (3=6.8 on 32 4 lattice indicating A/3 = 0.53± 
.02 (see Ref.5), we must stay at (3 < 6.87 ± 0.05 



Fig. 4(b) Ti n t a s a function of (3 at K=9. 
on a 32 4 lattice. 
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